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ABSTRACT 

Coherent  Anti-Stokes  Raman  Spectroscopy  (CARS)  is  a  new 
type  of  Raman  Spectroscopy.  The  phenomenon  is  associated  with 
the  nonlinear  conversion  of  two  laser  beams  in  a  medium  into  a 
third  collimated  beam  at  the  anti-Stokes  frequency  of  the 
medium.  This  analysis  evaluates  the  performance  capabilities 
of  CARS  under  various  experimental  conditions  for  the  molecules 
O2/N2/H2  and  CO.  Laser  intensity  fluctuations  and  shot  noise 
are  introduced  as  noise  sources  to  predict  signal-to-noise 
ratios  (S/N) .  The  S/N  is  evaluated  as  a  function  of  the  par¬ 
tial  pressure  of  a  gas  to  measure  the  performance  of  CARS.  An 
analysis  of  the  uncertainty  in  measuring  the  rotational  tempera¬ 
ture  is  presented  for  diatomic  gases.  Inverse  bremsstrahlung 
in  plasmas  is  investigated  as  a  possible  process  by  which  the 
parameters  measured  by  CARS  might  be  altered.  An  experiment 
is  discussed  in  which  an  attempt  was  made  to  remove  the  back¬ 
ground  signal  inherent  in  most  CARS  experiments.  Partial 
interferometric  cancellation  of  two  CARS  signals  was  observed 
by  the  proper  positioning  of  two  sample  cells  in  an  otherwise 
standard  CARS  experiment. 
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I. 


INTRODUCTION 


Coherent  Anti-Stokes  Raman  Spectroscopy  (CARS)  is  a  method 
of  chemical  analysis  yielding  molecular  spectroscopic  informa¬ 
tion  by  utilizing  nonlinear  scattering  phenomena.  As  a  part 
of  the  growing  field  of  nonlinear  optics,  CARS  has  been  used 
to  measure  nonlinear  coefficients  in  solids,  liquids, 
and  gases  [1-14].  Much  of  today's  efforts  are  focused 
on  the  use  of  CARS  as  a  nonreactive  probe  in  environments  not 
conducive  to  other  techniques  for  spectral  analysis  [14] . 

The  intent  of  this  thesis  is  to  assess  the  usefulness  of 
CARS  as  an  analytical  tool  using  the  example  of  diatomic  gases . 
CARS  is  related  to  the  more  basic  phenomenon  of  Raman  Scatter¬ 
ing  (RS) .  Raman  scattering  involves  a  shift  in  the  frequency 
of  a  small  portion  of  the  photons  incident  on  a  Raman  active 
medium.  RS  has  been  utilized  extensively  to  obtain  molecular 
spectroscopic  information.  CARS  yields  the  same  type  informa¬ 
tion  as  RS  except  that  the  collimated  CARS  signal  is  much 
easier  to  detect. 

CARS  was  feasible  for  chemical  analysis  only  after  the  ad¬ 
vent  of  tunable  lasers  of  moderately  high  power.  The  initial 
work  with  CARS  was  performed  by  Maker  and  Terhune  [15]  in  1963, 
but  only  recently  has  the  field  begun  to  expand  rapidly,  as 
evidenced  by  the  large  number  of  published  articles.  To  per¬ 
form  CARS  in  gases  requires  lasers  of  much  higher  power  than  in 
liquids  and  solids.  Typically,  pulsed  lasers  on  the  order  of 
megawatts  are  used. 
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The  use  of  such  high  power  pulsed  lasers  has  an  inherent 
problem  in  that  the  power  output  from  pulse  to  pulse  is  not 
constant.  This  introduces  a  source  of  noise  that  hinders  the 
observation  of  weak  resonances.  In  addition,  in  a  mixture  of 
gases,  a  CARS  signal  from  other  than  the  gas  of  interest  is 
present  as  a  background  signal.  These  two  factors  represent 
at  present  a  source  of  interference  in  observing  weak  molecular 
resonances . 

As  an  analytical  tool  CARS  can  yield  information  about 
molecular  identity,  concentration  and  rotational  temperature. 
Numerical  estimates  for  the  performance  capabilities  of  CARS 
are  based  on  the  ability  to  calculate  the  third-order  nonlinear 
dielectric  susceptibility,  This  susceptibility  is  com¬ 

posed  of  a  resonant  term,  y_,  and  a  nonresonant  term,  Xmd * 

Xo  is  calculable  from  the  physical  equations  that  define  it 
but  Xjjj^  was  deduced  from  experimental  observations. 

One  source  of  observations  for  XxTr>  relied  on  results  from 

^NR 

an  experiment  involving  Stimulated  Raman  Scattering  (SRS) . 
Appendix  A  contains  a  development  of  the  equations  related  to 
SRS  and  shows  the  relationship  of  the  SRS  gain  equation  to 
CARS. 

A  particularly  difficult  environment  to  perform  convention¬ 
al  spectral  analysis  with  Raman  scattering  is  in  an  ionized 
gas  called  a  plasma.  CARS  is  being  used  in  plasmas  and  is 
producing  better  results  than  conventional  Raman  scattering 
techniques.  An  analysis  is  made  to  predict  the  extent  of  in¬ 
teraction  between  a  plasma  and  the  laser  used  for  a  CARS 
experiment  in  the  plasma.  In  particular,  the  inverse 
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bremsstrahlung  process  for  the  transfer  of  energy  is  ana¬ 
lyzed. 

An  experimental  attempt  was  made  to  improve  the  basic 
CARS  experiment  by  the  elimination  of  the  background  signal. 
A  description  of  the  experiment  and  the  results  are  included 
in  this  study. 
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II. 


NATURE  OF  THE  PROBLEM 


The  intent  of  this  study  was  to  theoretically  investigate 
the  analytical  capabilities  of  CARS.  The  only  area  investi¬ 
gated  was  concerned  with  diatomic  gases.  It  was  desirable  to 
estimate  how  small  a  concentration  of  gas  could  be  detected 
and  with  what  accuracy  the  rotational  temperature  could  be 
measured.  CARS  has  demonstrated  its  usefulness  as  a  tool  for 
spectral  analysis.  It  offers  orders  of  magnitude  greater 
sensitivity  than  Raman  scattering  for  low  pressure  gases.  [16] 

Typically,  CARS  experiments  are  performed  by  focusing  two 
pulsed  laser  beams  of  different  frequencies  into  a  sample  cell. 
The  beams  cross  at  an  angle  that  corresponds  to  the  momentum 
conservation  requirement.  For  gases  this  angle  is  sufficient¬ 
ly  small  that  the  experiment  may  be  performed  with  the  beams 
colinear.  One  of  the  lasers  is  tunable  so  that  the  difference 
in  the  frequency  of  the  two  may  be  changed.  When  the  differ¬ 
ence  in  the  frequency  becomes  close  to  an  appropriate  natural 
molecular  frequency,  the  CARS  signal  increases  markedly  in 
amplitude . 

The  intensity  of  the  CARS  signal  is  dependent  on  the  pro¬ 
duct  of  the  square  of  the  intensity  of  one  beam  and  the  first 
power  of  the  intensity  of  the  other.  Any  amplitude  instability 
in  the  beams  leads  to  amplitude  instability  in  the  CARS  signal. 
Also,  the  presence  of  an  undesirable  background  signal  masks 
very  small  signals  of  interest. 
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A  brief  sketch  of  several  publications  in  the  field  of 
CARS  has  been  included  in  this  section  of  the  study.  The 
basic  equations  related  to  CARS  have  also  been  included  to 
develop  the  framework  on  which  the  numerical  calculations 
were  based. 

A.  BACKGROUND 

Maker  and  Terhune  [15],  in  1963,  were  the  first  to  perform 

a  CARS  experiment.  Prior  to  the  development  of  good  tunable 

laser  sources  most  of  the  work  with  CARS  was  developed  in  order 

to  measure  nonlinear  coefficients  of  materials. 

Yablonovitch,  Bloembergen  and  Wynne  [17]  were  among  the 

first  to  use  a  tunable  laser  to  generate  a  continuously  tunable 

difference  frequency.  They  investigated  solid  state  materials 

such  as  InSb  and  GaAs .  Levenson,  Flytzanis  and  Bloembergen 

[1]  followed  closely  with  an  experiment  in  diamond  to  measure 

(3) 

the  variations  in  X  as  a  function  of  the  difference  fre¬ 

quency.  DeMartini  [7]  utilized  CARS  with  H2  gas  in  1972  and 
Taran  [4]  in  1973  reported  that  he  had  used  CARS  as  a  means  of 
determining  the  spatial  distribution  of  H2  gas  in  a  flame. 

Taran 's  recent  work  has  resulted  in  the  spectra  of  several 
gases,  temperature  measurements  and  spatial  distributions  of 
gases  in  flames  [14] . 

B.  NONLINEAR  EFFECTS  IN  ISOTROPIC  MEDIA 

Due  to  the  nonlinear  dielectric  property  of  materials  two 
laser  beams  may  be  mixed  in  a  medium  to  generate  a  third  co¬ 
herent  beam  such  that 
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(1) 


0)2=  2  o)j^  -  0)2 

where  o)^  is  the  anti-Stokes  frequency  generated  from  the 
pump  and  Stokes  frequency  and  0)2  respectively). 

The  polarization  vector  associated  with  a  material  may  be 
expressed  as  a  power  series: 


P  (O)) 


E  (o))  + 


X 


(2)  ^2 

(O))  (O)) 


+  X 


(3)  -3 

(0))  ^(O)) 


+ .  .  . 


(2) 


where  X^^^  is  the  dielectric  susceptibility  tensor  of  rank 

i+1  and  E  is  the  applied  electric  field.  The  second  order 
(2 ) 

term  x  is  equal  to  zero  for  isotropic  liquids  and  gases 

because  of  inversion  symmetry.  The  lowest  order  term  that  is 

(3) 

present  in  a  gas  is  therefore  x  /  the  third-order  nonlinear 
susceptibility.  There  are  several  physical  processes  that  in¬ 
volve  the  third  order  term,  e.g.  ,  2a)^  + 

2a)j^  -  (1)2 etc.  The  last  process,  2a)j^  “  ^2'  pro¬ 

cess  referred  to  as  CARS. 

By  assuming  all  fields  are  along  one  axis  the  components 
of  the  tensor,  Xt  and  the  electric  field  vector  may  be 
treated  as  scalars.  The  magnitude  of  the  electric  field  may 
be  expressed  as : 


i  (k .  z 


Ei(a))  =  2 


-  0).  t) 
1 


+  c. c. ] 


(3) 


where  is  the  amplitude,  k^  is  the  propagation  constant 

(equal  to  m^nj^/c)  ,  2  is  the  length  along  the  z  axis  and 

t  is  time  in  seconds.  For  the  CARS  process  the  term  E^  in 


equation  (2)  includes  SE^ 
associated  with  CARS. 


E2  [18].  The  other  terms  are  not 
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The  resulting  third  order  polarization,  P 


(3)  . 

IS  given 


by ; 


p(3)_l 

“8 


(3) 

3x 


2  * 
01  02  e 


i  [  (2k^-k2)  z-  (2w^-a)2)  t] 


+C.C. 


(4) 


The  factor  of  3  has  been  a  point  of  inconsistency  among  au- 

(3) 

thors  in  the  past.  The  notation  for  x  (~“  /W,  ,a)  /W, )  is 

d.  iD  O  d 

consistent  with  that  of  Bloembergen  [19]  for  the  process 


0) 

a 


+  w  + 


w 


(5) 


Solving  Maxwell's  equations  for  plane  waves  with  polariza- 

(3)  .  .  . 

tion  P  the  anti-Stokes  gain  equation  is : 


iiTto^ 

■JciTJ 


i {2k^ -k  -k.) z 
)e  12  3 


(6) 


2 

In  terras  of  intensity,  where  |  ,  after. integrating 

equation  (6)  over  a  length  L,  the  result  is 


I 


3 


/  4tt  CO3 


AkL/z 


(7) 


where  Ak  =  2k--k.--k  .  The  conversion  efficiency  is  defined 

P3 

as  and  if  Ak  =  0  the  conversion  efficiency  becomes  a 

2 

maxiimom.  Because  of  the  small  dispersion  of  light  in  gas, 
CARS  may  be  performed  with  ,  ^2  and  colinear.  For  a 

focused  beam  (assumed  to  be  a  cylinder  of  plane  waves  and 
Ak  =  0)  the  conversion  efficiency  becomes: 


13 


The  third  order  nonlinear  dielectric  susceptibility, 
is  composed  of  a  frequency  dependent  resonant  part  (Xj^)  ^ 

nearly  frequency  independent  part,  the  nonresonant  susceptibili- 
tyr 

X  “  ^NR' 


X„  can  be  expressed  [18]  as: 
R. 


Xt5 


2Nc 


R  oi  4 
Jnoi, 


dn 


0) 


V 


[uj 


V 


-  iy  (a)^-(ii2)  1 


(10) 


When  a)^-u)2 


0)^,  the  resonant  susceptibility  becomes 


X 


R 


2Nc^  da  i 
3*10.2“ 


(11) 


where  N  is  the  number  density  of  molecules, 

Raman  cross  section  for  scattering,  A  is  the  difference  in 
probability  of  finding  a  molecule  in  the  ground  and  first 
vibrational  states,  and  y  is  the  line  width,  full  width  at 
half-maximum  intensity. 
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III.  NUMERICAL  ESTIMATION  OF 

(3) 

As  defined  in  Section  II. B.,  x  is  composed  of  two 

parts,  the  resonant  and  nonresonant  susceptibility.  The 
resonant  part  is  associated  with  the  various  Raman  active 
frequencies  that  characterize  a  molecule.  The  equation  for 
X„  has  been  derived  [18]  from  considerations  of  a  driven 
harmonic  oscillator  and  quantum  mechanics.  Xjgj^  is  associated 
with  the  electronic  structure  of  the  molecule  and  is  nearly 
frequency  independent. 

Equation  (10)  defines  Xc  i^  terms  of  an  entire  vibra- 

X\ 

tional  band,  namely  the  ground  state.  To  estimate  y-d 

R. 

each  rotational  level  requires  an  analysis  of  each  of  the 

\ 

terms  in  equation  (10)  that  are  functions  of  the  rotational 
quantum  number  J.  As  will  be  shown  in  Section  IV,  the  signal- 
to-noise  ratio  for  a  CARS  experiment  is  a  function  of  the  peak 
value  of  Xj^*  Being  able  to  predict  signal-to-noise  ratios  is 
important  when  deciding  which  analytical  technique  to  employ 
for  a  certain  situation.  Due  to  the  complexity  of  the  calcu¬ 
lations  a  computer  program  was  utilized  to  predict  the  peak 
value  of  X  for  several  gases. 

K. 

A.  RESONANT  SUSCEPTIBILITY,  Xo 

R 

1.  Cross  Section 

The  Raman  cross  section  for  scattering  is  typically 
quoted  as  the  cross  section  for  an  entire  vibrational  band. 

For  computational  purposes  it  was  necessary  to  apportion  this 
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total  value  among  the  rotational  levels .  The  relative  popula¬ 
tion  of  a  rotational  level  was  used  as  the  apportioning  factor 
such  that 

[  ^  (2J+1)  exp[-BJ(J+l)/kT]  / 

'^J  Q  \  dfi 

where 

Q  =  Z  g_  (2J+1)  exp [-BJ (J+l)/kT] , 

J 

J  is  the  rotational  quantum  nximber,  g  is  the  nuclear  spin 

u 

weighting  factor  and  B  is  the  rotational  constant.  The  values 
used  for  ^  were  from  Chang  [20] .  His  measurements  of  cross 
section  were  the  results  of  experiments  using  a  ruby  laser  at 
5140A.  The  cross  section  has  a  fourth  power  frequency  depend- 

4 

ence  (w  )  thus  an  accurate  calculation  of  X|^  requires  that 
the  data  be  corrected  for  this  factor. 

2 .  Linewidths 

On  resonance  the  CARS  signal  is  proportional  to  the 
inverse  square  of  the  linewidth.  Two  general  extremes  of 
linewidth  were  used,  the  Doppler  broadened  and  collision 
broadened  cases.  The  Doppler  broadened  lines  are  narrower 
than  the  collision  (or  pressure)  broadened  lines. 

3 .  Rotational  Transition  Frequencies 

For  molecules  that  exhibit  Raman  active  vibrational 
and  rotational  frequencies  the  CARS  signal  increases  in  ampli¬ 
tude  as  the  difference  frequency  approaches  one  of  the  molecu¬ 
lar  resonance  frequencies.  The  rotational  frequency  for  a 
vibration- rotation  transition  of  a  diatomic  molecule  is  ex¬ 
pressed  as: 


) 


(12) 


Total 
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(13) 


0)  T  =  03  -  203  X  (v  +  1)  -  a  J  (J  +  1) 

v,J  e  e  e  e 

where  03^  is  the  harmonic  vibrational  angular  frequency , 
is  the  enharmonic  constant,  v  and  J  are  the  quantum  numbers 
of  vibration  and  rotation  respectively,  and  is  the  rota¬ 

tion-vibration  interaction  constant. 

4 .  Nuclear  Spin  Weighting  Factor 

Homonuclear  diatomic  gas  molecules  have  a  statistical 
correction  factor  for  the  population  of  rotational  levels  due 
to  the  spin  of  each  of  the  nuclei.  Alternating  lines  in  the 
rotational  spectra  have  intensity  ratios  1;  (I/I+l)  where  I 
is  the  spin  of  each  nucleus.  The  strong  line  of  0^/02  and 
N2  are  the  even  numbered  J  levels  while  H2 ' s  strong  lines 
are  associated  with  odd  numbered  J  levels.  TABLE  1  summar¬ 
izes  the  degeneracy  weighting  factor,  g 

u  • 


TABLE  1 


Gas 


NUCLEAR  SPIN  WEIGHTING  FACTORS  (g_) 

\J 


Nuclear 

Spin 


(I) 


9j 

Even  #  J  Levels 


9j 


Odd  #  J  Levels 


22- 


1/2 


in 


^2- 


1/2 


Note  that  CO  is  not  a  homonuclear  molecule  and  thus  has  no 
degeneracy  factor;  i.e. ,  g,  =  1  all  J. 

u 

5 .  Computer  Program 

To  calculate  Xd  ^  nearly  continuous  function  of 

K. 

the  difference  frequency  required  a  computer  program.  A 


17 


Fortran  program  was  written  (See  page  54)  to  calculate  the 
maximum  amplitude  of  Xj^  ^  function  of  the  difference 
frequency  according  to  the  following  equation: 


.2Nc  da  A 


R,  (o)^  0)2) 


dQ  Q 


y ; v»»T 

“v,J  • 


(14) 


Tables  2  and  3  give  the  results  of  the  computations. 


B.  NONRESONANT  SUSCEPTIBILITY 

'^NR 

In  1967  Rado  [21]  performed  a  series  of  experiments  to 
measure  the  magnitude  of  the  nonresonant  susceptibility, 
of  several  gases.  The  impetus  for  the  experiments  was  to  mea¬ 
sure  the  peak  resonant  anti-Stokes^  power  from  a  hydrogen  cell 
and  compare  that  to  the  nonresonant  anti-Stokes  power  generated 
in  a  sample  cell.  The  ratio  of  powers  would  be  equal  to  the 
square  of  the  ratio  of  the  susceptibilities; 


sample  ^ 


'X  (Sample) y 
NR  J 

x'V 

P 


(15) 


The  value  Rado  chose  to  use  for  Xp  from  an  earlier 

work  with  phonon  lifetimes  and  SRS  gain  coefficients.  The 
relationship  of  Xp  the  SRS  gain  coefficient  is  discussed 


in  Appendix  A.  His  calculated  value  of 
.3 


'^NR 


for  N2  is 


1.35x10  .  Based  on  Table  3,  Xt^  N_  at  300  K 

erg  ,  ^P  2 


-16  cm' 


is  2.4x10 

erg 

susceptibilities  is 


Thus ,  the  ratio  of  the  peak  to  nonresonant 
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TABLE  2 


Gas 


^2 


«2 


CO 


Linewidth  (cm 

.  0034 

.0076 

.0055 

.012 


.  036 


.07 


.  005 


.07 


Teinp(°K)  Xp(cm  /erg) 


300 


8.3x10 


-15 


1500 


2. 0  X  10 


-16 


300 


3.3x10 


-15 


1500 


1  J.  X  10 


-16 


300 


1. 0  X  10 


-14 


1500 


4.  3  X  10 


-16 


300 


3.3  X  10 


-15 


1500 


2.2  X  10 


-17 
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TABLE  3 


Gas 

Linewidth  (cm 

Temp  ( ®  K ) 

Xp (cm^/erg) 

.16 

300 

2.6 

K  10-1« 

O2 

.07 

1500 

2.5 

X  10-1’ 

.14 

300 

2.4 

X  10-1« 

^2 

.06 

1500 

2.6 

X  10-1’ 

.03 

300 

1.2 

X  lO"^"* 

«2 

.08 

1500 

4.0 

in-16 

X  10 

.15 

300 

2.6 

X  10-1« 

CO 

.011 

1500 

1.1 

X  10-1^ 
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178. 


An  independent  method  to  estimate  the  magnitude  of 
was  provided  by  the  spectra  of  published  by  Taran  [14] . 

For  the  unresolved  spectra  the  difference  in  frequency  of  the 
maximum  and  minimum  values  of  the  CARS  signal,  6a),  is  given 
by  the  relation 


6a)  =  yj  1  + 


Xp2  \2 


(16) 


4X 


Assuming  a  linewidth  of  1.8  cm 
the  calculated  result  is  that 


NR 

-1 


and  a  6a)  of  42  cm 


-1 


47. 


(17) 


By  estimating  the  level  of  the  background  signal  in  the 
N2  spectrum,  another  ratio  of  equal  to 


N. 


75. 


(18) 


Knowing  the  peak  value  of  Xt,  yields  an  estimate  of 

R  NR 

An  earlier  work  by  Taran  [4]  provides  another  means  of 
estimating  the  ratio  of  the  susceptibilities.  The  total 
pressure  of  a  H2,N2  mixture  was  maintained  at  one  atmosphere. 
The  anti-Stokes  intensity  generated  was  plotted  as  a  fiinction 
of  the  H2  concentration. 
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scale  linearly  with  pressure,  then  the 


If  Xr  and  Xnr 
total  susceptibility  of 


the  mixture  may  be  approximated  by 


x<^’  =  XrP  + 

where  P  is  the  partial  pressure  of  H2.  As  a  ratio  of  Xj^ 

to  Y  the  observed  anti-Stokes  power  is  approximated  by 
'^NR 


I 


3 


(20) 


When  the  log  of  intensity  is  plotted  against  the  log  of  P 

(See  Figure  1) ,  the  predicted  and  observed  curves  are  in  fair 

4 

agreement  for  a  ratio  of  10  :1.  That  is 


X 


(H^) 


(21) 


Based  on 
of  Xp  for 
estimated  to 
ratio  is 


the  analysis  of  the  H^/N^  plot  and  using  the  value 
H-  from  Table  3,  the  value  of  Xt.tt^  for  N_  is 

2  2  2 

be  1.2x10  Thus,  another  estimate  of  the 

erg 


Xp  2  4xl0~^^ 

P  =  — --_ip  200. 

^NR  1.2x10 

The  factor  of  four  between  the  high  and  low  estimates  of  the 
ratio  is  not  surprising  based  on  the  variety  of  sources  used 
to  estimate  them. 
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The  solid 


Figure  1 


Taran's  results  follow  dotted  line 

X 

line  is  the  predicted  curve  for  — 

X 

3 

two  curves  are  equal  from  10  to 


R 

NR 


=  lo"^.  The 
10^  ppm. 


H2  Concentration  (ppm) 


IV.  STATISTICAL  CONSIDERATIONS 


The  anti-Stokes  power  resulting  from  a  tightly  focused 
CARS  experiment  is  expressed  as  [18] : 


P3=K 


3X 


(3) 


2  T.  2„ 

P  P 
1  2 


(22) 


where 


K 


■fe) 


2  47T  03- 


(23) 


The  resonant  susceptibility  may  be  expressed  in  terms  of  its 

peak  value,  x  /  t»y; 

P 


XpY* 

^R  ~  2Ao3  -  iy 


(24) 


where  A03  =  03,  “03.,-03  . 

1  2  V 

For  weak  signals  in  which  x  absolute  value 


squared  of  x 


(3) 


is  given  by 


X 


(3) 


X-n  Y 


=  lx 


NR  2Ao3-iy 


(3) 


%R 


2xpYAo3 

2  2  ■ 
4 (Aw)  +Y 


(25) 


For  a 


The  magnitudes  of  x'^^  its  extrema  are  Xjgj^jlXp/2’ 

signal  strength  proportional  to  the  difference  in  the  peak  and 

background  susceptibility 


Xt.  2 


1 

2 


^•5“(Xnr  2  ^  ~  ^NR  ^NR^P 


^NR 


(26) 


Thus,  for  x-D  <<  Xmtd'  anti-Stokes  power  above  that  due 

IT  ^NR. 

to  Xj^R  alone,  is 


24 


(27) 


^  ^NR  ^1  ^2'^T 


where  rirj,  is  an  efficiency  coefficient  to  take  into  account 
the  optical  losses,  detector  inefficiencies  and  observed  con¬ 
version  efficiency  factor. 

Typically,  the  largest  source  of  noise  when  performing 
CARS  on  samples  at  one  atmosphere  pressure  is  associated  with 
the  power  fluctuations  of  the  pump  and  Stokes  beams.  If  the 
power  of  the  laser  source,  ,  has  a  relative  standard  devia- 
tion  of  e  ,  the  standard  deviation  of  the  power  is  e  P^ .  The 

P  P  Li 

fluctuation  of  the  background  signal  power  P^,  due  to 
is  thus : 


'’n  “  ^  ’‘nr 


(28) 


assuming  each  beam  has  the  same  relative  standard  deviation  e^. 

For  a  low  intensity  CARS  signal  in  the  absence  of  back¬ 
ground  one  finds  the  major  source  of  noise  to  be  the  shot  noise 
of  the  signal.  For  an  N  photon  signal  the  standard  deviation 
is  /W. 

For  a  pulsed  signal  the  noise  power  due  to  shot  noise  is 


2  ^3^^ 


N 


(29) 


The  signal-to-noise  ratio  for  the  background  limited  case 
given  by  equations  (27)  and  (28)  becomes 


S/N 


P _ 

%R 


(30) 


whereas  the  signal-to-noise  ratio  for  a  situation  in  which 
there  is  only  a  very  weak  signal  is  given  by  equations  (27) 
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and  (29)  to  be 


1  1 

S/N  (K  Xp^  (T/hv)^.  (31) 

As  an  application  of  the  signal-to-noise  problem  consider 
the  measurement  of  Xp  from  a  spectral  envelope.  The  uncer¬ 
tainty  of  Xp,  Op,  is  expressed  as;  [22] 


a  ^  ^  a. ^ 

P  1 


(32) 


where  is  the  standard  deviation  in  determining  tp  from 

th 

the  i  pulse  of  a  pulsed  laser.  For  CARS  with  a  background 

limited  situation  and  uncertainty  in  the  signal 

amplitude  is  approximately  a  constant,  associated  with 

the  laser  power  fluctuations. 

3Xp 

Thus  =  a_  ( - )  where  a  is  basically  the 


m 


hi 


m 


(w.  J,T) 
1  / 


noise  power  given  by  equation  (28): 


a  =  3  e  K  xlry 
m  p  '^NR  1  2 


(33) 


I,  T  is  the  intensity  of  the  spectral  envelope  at  co.  , 

((Jd  ,  U  f  1  )  1 

1  / 

due  to  the  quantum  rotational  number  J,  at  temperature  T. 
The  result  of  this  analysis  yields  [16] : 


a 


a  = 


m 


/ 2Aaj>  2 
'  TTY 


(34) 


or 

1 

_  1  .2Aa).  2 

^p  S/N  ^  TTY 


(35) 
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In  general  any  standard  deviation  is  dependent  on  the  nuniber 
of  times,  N,  the  measurement  is  taken  such  that 


a 

P 


1 


a  — 


The  conclusion  is  that  for  a  given  S/N  the  term  is  the 

effective  number  of  pulses  required  at  intervals  of  Ao)  to  pro 
duce  an  uncertainty,  a^,  in  the  measurement  of  Xp- 

The  uncertainty  in  determining  the  rotational  temperature 
is  a  product  of  the  nearly  constant  background  uncertainty, 
and  the  small  change  of  intensity  with  temperature.  A 
convenient  definition  for  a  is: 


m 


a 


m 


I  (5, 
max  I 


(36) 


max 


where  6^  is  the  relative  uncertainty  in  the  amplitude  of 
■^max 

the  line  of  maximiim  intensity.  The  uncertainty  in  the  tempera 
ture  then  becomes  [16] : 


1 


(37) 


The  intensity  of  a  line  is  proportional  to  the  Boltzmann 
population  distribution.  The  change  of  intensity  with  tempera 
ture  has  been  shown  [16]  to  be; 


^^J,T  ^  [BJ(J  +  1)  -  kT] 

3T  J,T  ^^2 

The  resultant  standard  deviation  in  temperature  is 
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It  has  been  estimated  [16]  that  for  an  uncertainty  in  the 
maximum  intensity  equal  to  2%  the  temperature  of  CO  at  1500®K 
may  be  determined  to  within  10 °K. 

This  uncertainty  in  temperature  is  for  an  unrealistic 
situation  in  that  each  rotational  line  is  completely  re¬ 
solved.  A'  more  realistic  situation  is  the  case  of 
a  continuous  envelope  as  given  in  equation  (14) .  A  computer 
program  (see  page  59)  was  written  to  calculate  the  uncertainty 
in  a  temperature  measurement,  based  on  equation  (37). 


,  31(0).  J,T)  2 

y  ( _ ij: _ ) 

^  2  i  ^  3T  ^ 

m 


o  =  Ie  =  KY_e 
m  p  p  ''•P  p 


,  3(K  J,T) 


For  the  same  conditions  used  to  calculate  the  previous 
uncertainty  in  temperature  the  results  of  the  envelope  calcu¬ 
lation  yields  an  uncertainty  of  15 °K.  This  value  is  of  the 
same  order  of  magnitude  as  the  statistical  prediction  of  10°K, 
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estimated  by  using  equation  (39) .  It  is  to  be  expected  that 
the  uncertainty  would  be  higher  for  measurements  on  the 
spectral  envelope.  Verification  of  the  statistical  approach 
means  that  the  much  simpler  method  for  predicting  the  uncer¬ 
tainty  is  a  good  approximation. 
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■  V.  CARS  IN'  A  PLASMA 


Low  pressure  ionized  gases,  or  plasmas,  represent  a  class 
of  materials  that  are  receiving  very  close  investigation.  The 
population  of  the  various  energy  levels  reveals  much  about  the 
plasma.  Determining  the  population  of  levels  is  a  very  diffi¬ 
cult  task  with  conventional  Raman  spectroscopy.  CARS  has  been 
utilized  to  obtain  spectra  of  several  plasmas  at  the  Naval  Re¬ 
search  Laboratory,  Washington,  D.  C.  The  experience  there 
indicates  that  the  technique  is  not  difficult  and  produces 
useful  information. 

Of  importance  when  analyzing  CARS  in  a  plasma  is  how  much 
the  high  power  laser  beams  may_  be  expected  to  perturb  the 
system  under  investigation.  The  major  source  of  interaction 
between  the  plasma  and  laser  beams  appears  to  be  from  the  in¬ 
verse  bremsstrahlung  process  [23]  .  The  free  electrons  in  the 
plasma  will  gain  kinetic  energy  by  this  process. 

If  the  laser  frequency  is  above  the  plasma  frequency  the 
light  will  propagate  through  but  there  will  be  some  attenuation 
or  absorption.  The  plasma  frequency  is  given  by  [24] 


1 


-3 

where  ri^  is  the  number  density  of  electrons  (cm  )  ,  e 
is  the  charge  of  an  electron  (esu)  and  is  the  mass  of  an 

electron  (grams) . 

An  estimate  of  the  change  in  the  electron  kinetic  energy 
may  be  made  by  assuming  there  is  an  absorption  coefficient,  a. 
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for  the  plasma.  The  energy  absorbed  per  unit  volume  of  the 
plasma  is 


"ABS 


T-  ,,  -ox,  ^ 
I  (1  -  e  )  At A 
o _ 

Ax 


P^Ata 

o 


(41) 


where  I  is  the  incident  intensity,  P  is  the  incident 
o  .  o 

power.  At  is  the  pulse  duration,  A  is  the  beam  cross  sec¬ 
tion  and  X  is  the  path  length  of  the  medium.  If  all  the 
absorbed  energy  is  transferred  to  the  kinetic  energy  of  the 
electrons  then  the  energy  absorbed  per  electron  is 


aP  At 

E  —  ^ 

ABS  An 

e 

As  a  fraction  of  the  initial  energy,  3/2kT^, 


^ABS 

3/2kT^ 


aP  At 
o 


3/2An^kT^ 


(42) 


(43) 


where  is  the  electron  temperature.  The 'absorption  co¬ 

efficient,  a,  for  inverse  bremsstrahlung  has  been  reported 
as  [25] 

1  ^2  6 

®  ^  ^3/2kT^  r  3/2  3  “  exp  (-j^)  ]  (44) 

h  c  M  '  V 
e 

where  n^  is  the  ion  number  density,  z  is  the  charge  of  the 
ions  and  v  is  the  frequency  of  the  incident  light.  Evalua¬ 
tion  of  equation  (43)  with  a  defined  by  equation  (44)  results 
in  a  relative  change  of  temperature  equal  to 

3T^  1.3xl0^^n^z^P^At[l  -  exp(-hv/kT)] 


T 


3 


3/2 


(45) 


This  agrees  favorably  with  an  equation  for  the  change  of 
electron  temperature  given  by  [23] 


25  2 

AkT  1.4x10  n • z  [1  -  exp (-hv/kT) ] r  P  dt 

_ e  _  _ 1 _ 0 _ -i-i 


kT 


3,10,  x3/2 
V  (kT  )  ' 
e 


(46) 


where  electron  volts  and  is  the  power  of 

the  laser  in  watts. 

The  equipartion  time,  t^,  is  the  average  time  required 
for  the  electrons  to  transfer  their  excess  kinetic  energy  to 
the  ions  and  thus  reach  an  equilibrium  temperature.  The  equa¬ 


tion  for  t^  is  given  by  [25] 


t  = 
e 


25.2  A(T  ) 
e 


3/2 


(47) 


where  A  is  the  atomic  weight  of  the  ions ,  T^  is  the  electron 
temperature  in  degrees  Kelvin,  and  is  in  cm  If  the 

duration  of  the  laser  pulse  is  shorter  than  t^  then  any 
energy  the  electrons  gained  by  inverse  bremsstrahlung  will 
not  be  transferred  to  the  ions  and  there  will  be  very  little 
perturbation  of  the  molecules  under  investigation. 

Consider  the  following  as  a  typical  example: 

4 

plasma;  Deuterium  at  50  Torr,Te  =10  K 

11  -3 

7e  =  10  cm  ,  20  cm  long 
•  -4  2 

laser;  5500A,  5MW,  10  cm  beam  area,  20  nsec 


duration 

The  plasma  frequency, 

The  equipartion  time ,  t^ , 


equals  1.8  x  10^^  sec  ^  and  • 

L  p 

equals  20  micro  seconds.  The 


32 


relative  change  in  electron  temperature  given  by  equation  (46) 
is  .0001%.  Little  perturbation  of  the  plasma  is  expected  from 
these  conditions. 
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VI.  BACKGROUND  SUPPRESSION  EXPERIMENT 


As  discussed  earlier,  the  nonresonant  susceptibility,  Xmd' 
is  responsible  for  the  background  signal.  The  power  amplitude 
fluctuations  of  the  lasers  produce  noise  on  this  signal.  The 
resonant  signal  above  the  background  also  has  noise  on  it.  It 
would  be  much  easier  to  detect  the  resonant  signal  if  there 
were  no  noise  on  the  background  or,  better  yet,  no  background 
at  all.  The  desire  for  reduction  of  the  background  signal  led 
to  an  idea  [26]  for  an  experiment  to  determine  if  it  was  possi¬ 
ble  to  completely  suppress  the  background.  In  addition  to 
yielding  much  better  signal  quality  the  elimination  of  the 
background  would  limit  the  detectibility  of  CARS  to  the  limit 
of  conversion  efficiency. 

Typically,  the  diluent  of  a  sample  provides  the  largest 
contribution  to  The  basic  idea  for  suppression  of  the 

background  signal  was  to  perform  two  CARS  experiments  in  series . 
In  the  first  experiment  the  generated  anti-Stokes  beam,  co^/ 
would  be  due  to  only  ^  diluent)  and  in  the  other 

I 

experiment  the  anti-Stokes  beam,  would  be  due  to 

Xj^  (the  sample  plus  the  diluent)  .  If  the  to^  beam,  from  the 

first  cell,  were  then  forced  out  of  phase  by  Ntt  radians 

I 

(where  N  is  an  odd  integer)  relative  to  the  second  beam, 
and  superimposed  on  the  second  then  the  resultant  signal  should 
be  only  that  due  to  Xo  of  sample. 

Molectron  Corporation  of  Sunnyvale,  California,  through 
Dr.  G.  Klauminzer  provided  equipment,  facilities  and  assistance 

for  conducting  the  experiments . 
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Figure  2.  Arrangement  of  Equipment  for  Background  Suppression 
Experiment. 

Figure  2  shows  the  basic  arrangement  found  to  be  most 
convenient.  A  pulsed  nitrogen  laser  was  used  to  drive  two 
tunable  dye  lasers  in  an  arrangement  reported  by  Dr. 

Klauminzer  in  Reference  12. 

Benzene  was  chosen  as  the  first  material  to  work  with 
since  it  was  known  that  the  CARS  signal  was  very  strong  and 
visible.  For  benzene  the  values  for  W2  and  were 

480.,  504.  and  458.2  nm,  respectively. 

Both  cells  were  filled  with  benzene  in  anticipation  of 

f 

demonstrating  complete  cancellation  of  with  co^-  It  was 
known  that  a  translation  of  cell  one  along  the  focal  axis 
would  change  the  intensity  of  by  changing  the  volume  in 

the  focal  region  but  it  was  also  discovered  that  the  same 

I 

translation  changed  the  phase  of  relative  to  To 

f 

completely  cancel  it  was  known  that  several  conditions 

had  to  be  met.  The  beam  had  to  be  ir  radians  out  of  phase,  the 
intensities  had  to  be  equal  and  the  beams  had  to  be  coinci¬ 
dent.  A  plot  of  the  resultant  intensity  versus  the  translation 
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of  cell  one  revealed  that  there  was  constructive  and  destruc¬ 
tive  interference  of  the  two  beams.  The  destructive  interfer¬ 
ence  was  not  complete  though.  A  visual  observation  of  the 

I 

co^  +  beam  revealed  interference  fringes .  By  steering  the 

I 

beam  off  of  the  beam  the  fringe  pattern  could  be 

changed.  It  was  suspected  that  chromatic  aberration  in  the 
focal  lens  was  responsible.  Calculation  of  the  effect  of 
chromatic  aberration  indicated  that  the  fringe  pattern  should 
be  approximately  that  observed.  The  use  of  an  achromat  for 
the  focal  lens  reduced  the  number  of  fringes  visible  but  it 
was  found  that  a  slight  rotation  of  the  lens  changed  the 
pattern  considerably. 

The  best  measured  ratio  of  minimum  to  maximum  intensity 
as  a  function  of  cell  translation  was  observed  to  be  10:1. 
Since  it  was  the  electric  fields  that  were  interfering,  the 
10:1  ratio  amounts  to  about  a  60%  cancellation  factor.  The 
cancellation  factor  was  arrived  at  by  considering  the  incom¬ 
plete  phase  cancellation  of  the  electric  fields.  The  observed 
intensity  is  proportional  to 


+  E, 


i0 


e 


2 


2 


(48) 


If  the  amplitudes  and  are  equal,  the  resulting 

intensity  is  proportional  to 

P  ie  i0„  -i0,  -i0„ 

E  (e  -^  +  e  (e  -^  +  e  (49) 


or 


a 


2 

E  (2 


+  2  cos  (0, 


-  02)) 


(50) 
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results 


Complete  constructive  interference,  when  0^^  =  Q^, 
in  an  intensity  four  times  that  of  one  beam,  and  complete  de¬ 
structive  interference  would  result  in  no  intensity. 

It  is  not  surprising  that  on  the  average  complete  inter¬ 
ference  does  not  take  place.  For  a  small  difference  in  0^ 

2 

and  02  the  observed  magnitude  is  slightly  less  than  4E  or 
slightly  more  than  zero. 

An  observed  9:1  intensity  ratio  is  equivalent  to  a  3.6  to 

2 

.4  ratio  on  the  normalized  scale,  E  =1.  Thus,  the  inter¬ 
ference  produces  a  signal  equal  to  40%  of  what  one  beam  would 
produce,  therefore,  a  cancellation  factor  of  60%. 

Next,  a  three  molar  sodium  benzoate  solution  was  used  in 
the  cells  in  place  of  benzene.  Again  there  was  incomplete 
cancellation.  The  best  observed  ratio  of  intensities  was  5:1. 

It  had  been  predicted  that  sweeping  the  frequency  of  0^2 
record  a  spectriom  would  not  alter  the  phase  difference  of 

I 

and  The  first  observations  made  while  sweeping  0)2  re- 

I 

vealed  that  and  went  in  and  out  of  phase  over  a 

period  of  about  10  cm  By  changing  the  optical  path  length 

of  0)^  or  and  0)2  by  use  of  microscope  slides  the 

period  could  be  increased  to  about  60  cm 

The  observed  increase  in  the  period  of  the  oscillation 
verified  that  the  optical  path  length  could  be  controlled  to 

I 

shift  0)^  and  in  and  out  of  phase.  With  the  appropri¬ 

ate  optical  path  length  the  period  of  the  oscillation  could 
be  adjusted  to  infinity. 

A  spectrum  was  recorded  for  the  three  molar  sodium  benzoate, 
two  cell,  series  cancellation  experiment.  Of  interest  from  the 
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spectra  is  the  interference  of  the  background  with  respect  to 
the  noise  levels  observed.  Figure  3  is  a  portion  of  that 
spectrum  depicting  the  observed  resultant  signal.  When  the 
interference  was  constructive  the  noise  was  much  higher  than 
for  the  destructive  interference. 

As  stated  earlier,  the  reason  for  the  experiment  was  to 
increase  the  sensitivity  of  CARS.  An  analysis  of  the  signal- 
to-noise  ratio  for  a  cancellation  experiment  indicates  that 
the  cancellation  has  to  be  much  more  effective  than  was  ob¬ 
served  in  this  experiment  to  gain  an  order  of  magnitude  in¬ 
crease  in  the  signal- to-noise  ratio. 

2 

Assxime  b  is  the  fraction  of  background  intensity  remain¬ 
ing  after  cancellation.  Then,  the  S/N  ratio  is  given  by 


S/N  = 


t 


similar  to  the  development  of  equation  (30) .  An  increase  in 

2 

S/N  of  10  requires  b  =  .01,  or  a  99%  cancellation  factor. 

A  400:1  intensity  ratio  would  correspond  to  such  a  cancella¬ 
tion  factor. 

The  experiment  as  performed  was  very  sensitive  to  the 
alignment  and  quality  of  the  optics  involved.  One  of  the 
major  contributions  to  the  incomplete  cancellation  is  probably 
the  diffraction  limited  ability  of  the  center  lens  to  perfect- 

I 

ly  image  the  initial  tu^  beam  onto  Other  macroscopic 

factors  such  as  turbulence  and  inhomogeneity  of  the  samples 
are  other  obstacles  to  complete  cancellation. 
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background  signal 


VII.  RESULTS  AND  DISCUSSION 


The  basis  for  determining  the  areas  of  application  for  a 
new  analytical  tool,  such  as  CARS,  must  depend  on  recognizing 
the  limits  beyond  which  the  technique  cannot  be  expected  to 
yield  the  desired  information. 

The  results  of  this  analysis  are  presented  as  lower  limits 
of  concentration  that  are  expected  to  be  detectable  by  CARS. 
For  a  mixture  of  gases,  the  detectability  of  one  gas  depends 
on  the  ratio  of  Xo  Xmt,/  the  same  as  Taran's  plot  for 

H2  and  N2.  Plots  of  the  ratios  for  other  gases  with  N2  as 
a  diluent  are  presented  in  Figure  4.  The  minimxim  detectable 
concentration  is  dependent  on  the  required  signal-to-noise 
ratio  for  a  given  experimental  condition.  The  lower  limit  of 
detection  for  a  pure  gas  at  reduced  pressures  would  be  limited 
to  the  conversion  efficiency  if  not  for  the  photon  statistics 
involved.  From  the  signal-to-noise  analysis,  estimates  of  the 
detectable  lower  pressures  are  summarized  in  Figures  5  and  6. 

Based  on  the  results  of  Section  IV  the  rotational  tempera¬ 
ture  of  a  gas  is  expected  to  be  measurable  to  within  1%.  That 
assumes  the  amplitude  fluctuations  of  the  CARS  signal  observed 
remain  less  than  10%. 

In  a  plasma,  the  inverse  bremsstrahlung  process  is  ex¬ 
pected  to  be  a  relatively  insignificant  perturbation  to  the 
parameters  determined  by  a  CARS  experiment. 

The  experimental  attempt  at  cancellation  of  the  background 
confirmed  the  idea  that  interferometric  cancellation  is  a 
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possible  technique  for  improving  CARS  when  the  background  be¬ 
comes  a  problem.  The  observed  results  were  not  of  the  magni¬ 
tude  hoped  for  but  did  serve  to  verify  the  concept  and  to 
recognize  a  possible  source  of  limiting  performance. 

Turbulence  in  the  samples  and  inhomogeneity  may  be  part 
of  the  reason  that  incomplete  cancellation  occurs.  The  major 
contribution  to  the  lack  of  complete  cancellation  is  probably 
the  result  of  diffraction  limited  optics.  The  first  beam 
cannot  be  refocused  through  a  lens  and  have  precisely  the  same 
geometry  that  it  had  originally. 
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Relative  Anti-Stokes  Intensity 


I 


Figure  4 


Concentration  [ppm]  of  Sample 


Log-lo^  plots  of  N_  gas  mixtures  at  1 
total  pressure.  at  300  K.  Upper  CO 

300  K,  Lower  coplot  at  1500  K.  of 

to  be  1.2  X  lO"^®  at  300  K,  and 

-IS  3 

.24  X  10  at  1500  K. 
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atmosphere 
plot  at 
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Figure  5. 


Signal-to-noise  ratios  for  gas  mixtures  at  300  K 
with  a  total  pressure  of  1  atmosphere.  Background 
pressure  and  signal  due  to  N2.  The  area  boundar¬ 
ies  are  defined  by  assuming  values  of  the  un¬ 


certainty  in  the  pulsed  laser  intensity, 

-2  -1 
10  to  10  . 


£  ,  from 
P 
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Signal-to-Noise  Ratio 


PRESSURE  (ATMOSPHERES) 


Figure  6.  Signal~to~noise  ratios  for  Doppler  broadened 

gas  at  300®K  as  a  function  of  pressure.  Noise 
due  to  photon  statistics.  The  assumed  parameters 
are 

=  .01,  P^  =  5MtV,  P2  =  .5MW,  T  =  20  nsec, 

X  =  5500  A. 
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VIII.  CONCLUSIONS 


This  analysis  presents  a  model  by  which  the  limits  of  the 
information  available  from  CARS  may  be  predicted.  At  STP 
typical  diatomic  gases  should  have  a  signal-to-noise  ratio  of 

-3 

1  at  partial  pressures  of  about  10  atmosphere.  Hydrogen,  be 

cause  of  its  narrow  linewidth  and  population  distribution, 

“5 

should  have  a  signal-to-noise  ratio  of  1  at  about  10  atmos¬ 
phere.  Pure  gases  at  low  pressure  where  the  linewidth  is 
determined  primarily  by  Doppler  broadening,  may  be  detected 
at  pressures  several  orders  of  magnitude  lower  than  this. 

The  ultimate  detectability  depends  upon  laser  power  and  satura 
tion  effects. 

Simple  relationships  have  been  derived  which  allow  the 
prediction  of  the  uncertainty  in  temperature  of  a  gas  for  an 
assumed  signal-to-noise  ratio. 

In  a  plasma  the  inverse  bremsstrahlung  process  is  found 
to  transfer  insufficient  energy  to  the  plasma  to  significant¬ 
ly  modify  the  parameters  measured. 

The  experimental  attempt  at  background  reduction  by  inter¬ 
ferometric  cancellation  verified  that  the  principles  were 
correct.  Complete  cancellation,  however,  was  not  observed. 
This  is  thought  to  be  due  to  the  diffraction  limitations  of 
optical  lenses. 

One  area  of  additional  research  and  experimentation  that 
could  further  extend  the  applications  of  CARS  is  that  of  the 
background  suppression.  Another  area  of  improvement  needs  to 
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be  that  of  the  fluctuating  intensities  of  the  high  power 
pulsed  lasers  that  are  available. 
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APPENDIX  A 


STIMULATED  RAMAN  SCATTERING 


A.  DEVELOPMENT 

The  interaction  of  two  electromagnetic  fields  in  a  material 
is  such  that  different  combinations  of  the  two  may  produce 
additional  fields.  The  two  additional  fields  of  interest  here 
are  the  Stokes  and  anti-Stokes  frequencies. 

Raman  scattering  has  been  observed  extensively  and  found 
to  produce  the  Stokes  and  anti-Stokes  frequencies  at  plus  and 
minus  the  vibrational  frequency,  of  the  molecule  from 

the  incident  light,  03^. 


"'stokes  -  "'l  ■  "V 


"'anti-Stokes  "'l 


+  CO. 


V 


(Al) 

(A2) 


Stimulated  Raman  scattering  predicts  a  gain  in  the  strength 
of  a  field  as  it  traverses  a  media.  The  media  of  interest  in 
this  study  are  gases.  The  gain  will  be  related  to  the  non¬ 
linear  susceptibility  of  those  gases. 

Beginning  with  Maxwell ' s  equations ,  which  contain  various 
assumptions  implied  by  their  form,  the  plane  wave  equation  is 
derived  from  the  curl  of  the  curl  of  the  electric  field,  E. 
Assuming  a  solution  to  the  wave  equation  of  the  form: 


,  -i (kz 

E  .  i  (£  e 


cot) 

+  c. c. ) 


(A3) 


where  e  is  the  amplitude ,  k  is  the  propagation  constant 
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and  w  is  the  frequency.  For  simplicity,  all  fields  have 
been  reduced  to  one  dimension  along  the  z  axis.  From  the 
nonlinear  polarization  term  the  various  possible  frequencies 
are  defined  and  their  field  properties  (gain  or  attenuation) 
are  explained  by  computing  the  solution  to  the  wave  equation. 
For  completeness  both  MKSA  and  CGS  results  are  presented. 

The  development  in  MKSA  units  begins  with  Maxwell's  basic 
and  constituent  equations. 


V  X  H  =  J  +  D 


D  = 


e  E  +  P 
o 


V  X  E  =  -B 

V  •  E  =  0 


J  =  aE 


P  =  X  E  +  P 


NL 


B  =  0 


B  = 


The  curl  of  H  becomes 


V  X  H  =  aE  +  e^E  +  P  =  aE  +  eE  + 

Then  the  wave  equation,  defined  from  VX(VxE)  is 


-V 


dt 


(V  •  B) 


d 

"’^o  dt 


(V  X  H) 


or 


E  + 


V^E 

^o" 


£  NL* 


(A4) 


Using  the  assumed  form  of  the  solution,  equation  (A3) ,  and 

2  2 

the  assumptions  that  a  =  0  and  3  e/3z  is  ve3ry  small, 
then  equation  (A4)  becomes 
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i (kz  -  wt) 


0)  , 
(ee 


2  2 

.  „  ^  V  ,  k  c  , _ i(kz  -  a)t)  ^  . 

+  c.c.)  +  - y-  (ee  '  +  c.c.) 

2n 


kc' 


o  Z 


.c.j  = 


-  i  i 

e  NL' 


(A5) 


The  usual  definitions: 


2  1 
c  = 


u  e 
o  o 


r)  =  /e/e 

o 

k  =  0)  ri/c 

have  also  been  used. 

For  the  two  interactive  fields  and  the  nonlinear 

polarization,  P  is  expressed  as:  [18] 

NLj 


P  =  ^(3)„  2 

^NL  ^  o^  ^1  ^2* 


(A6) 


Three  frequencies  are  possible  from  the  E2  combination; 

0)2,  2a3^-aj2 / 210^+0)2 •  If  ^2  Stokes  field  with  frequency 

0)2  (assuming  aj^<  ^2)  then  a  resultant  field,  E^,  of  fre¬ 
quency  210^-032  is  the  anti-Stokes  field.  For  the  anti-Stokes 
calculations  is  given  by: 

o  *  -i  [  (2k, -k_)  z- (2a),-u3«)  t] 

^NL  ^  ^^o^  ^8‘^1^1^2  c.c.)].  (A7) 


The 


resultant  field  may  be  E2  for  the  case  when 


3  ,1 


'’nL  '  <3=1^1  ^2' 


-i  (k2Z  -  a32t) 


+  c.c. ] 


(A8) 


49 


1 .  Stokes  Gain  Equation 


To  solve  equation  (A5)  for  the  Stokes  condition,  equa¬ 
tion  (AS) ,  there  are  two  possible  alternatives  that  lead  to  a 
gain  equation.  One  way  is  to  assume  that  9 £2/22  is  very 
small  and  k2  is  complex  such  that 

k2  =  k'  +  ik". 

The  imaginary  part  of  k2,  k” ,  is  then  evaluated  as: 

U)  3x"lEi°l^ 

V »»  —  ^ _ i _ 

^  2cn 

where  x"  is  the  imaginary  part  of  is  the 

magnitude  of  at  z  =  0. 

The  observed  intensity  of  E2  is  proportional  to  the 
2 

time  average  of  IE2I  or 

where 


K  T.  *  1/  ik'2  -k"z  ^  *  -k'z  -k”z,  ,  ^  , 

2^2^2  “  8  ^2®  ®  ^2  ®  ®  (c.c.) 


Thus,  1^2 1  proportional  to  e  .  x”  is  a  negative 

— 2k*'z 

quantity  [27]  and  therefore  e  is  an  exponentially 

growing  function  with  distance. 

The  other  possible  way  to  solve  equation  (A5)  is  to 


let 


and  then  solve  for 


9e2/9z*  The  result  is  that 
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2 

3(X'  +  ix")  |Ei°| 


1 

2^2 


] 


and  again  for  x"  negative  there  is  gain  for  increasing  z. 

Both  approaches  result  in  an  intensity  gain  coeffici¬ 
ent,  g,  such  that 


I 


2 


I 


o  gz 
2  ® 


where  is  the  intensity  of  at  z 

that 


g  = 


|3X”I 


2 


0 .  The  result  is 


{A9) 


For  cgs  units  the  gain  is 


g 


477032 


n2c 


(AlO) 


2 .  Anti-Stokes  Gain  Equation 

Solving  equation  (A5)  for  the  anti-Stokes  condition 
(A?)  results  in  the  gain  equation  for  the  anti-Stokes  field. 
The  only  way  to  solve  equation  (A5)  is  to  assume 


The  result  is  that 


3^3  _  ^“3 

3z  n  3C 


3(X’  +  ix”) 


1  *  iAkz 

8^1^1^2  ® 


(All) 


where  Ak  =  2k^  -  k2  -  k^. 

Again  if  x"  is  negative,  there  is  gain.  In  cgs  units 
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(A12) 


ae. 


3e . 


3z  (cgs) 


=  4tt  (- 


3z  'MKSA. 


B.  CORRELATION  OF  VARIOUS  ARTICLES 

The  only  published  values  for  nonresonant  susceptibilities 
of  gases  are  those  given  by  Rado  [21].  It  was  desirable  to 
check  for  consistency  in  the  equations  used  by  Rado  and  the 
ones  used  to  calculate  resonant  susceptibilities. 

The  basis  for  Rado's  experimental  values  of  '’^he 

resonant  value  of  ^2  other  experimental 

data  [23] .  The  gain  coefficient  for  Stokes  Stimulated  Raman 
Scattering  provided  Rado  with  the  needed  value  of  Xn 
at  pressures  greater  than  10  atmospheres.  The  gain  coeffici¬ 
ent  was  given  in  [28]  in  cgs  units  as 


k [Stt^L  a  Na^]  T 
g  _  _ _ _  r 

I,  e  cn  "K 


X  10  (cm/W) 


{A13) 


where  L  is  the  local  field  correction  factor,  N  is  the 
number  density  of  atoms,  A  is  the  population  ratio,  is 

the  phonon  lifetime,  a  is  the  polarization  coefficient  and 
Ij^  is  the  incident  laser  intensity. 

An  independent  article  by  Maier  [29]  used  the  following 
definition  for  the  gain  coefficient 


g  _  8itN  (Jpg 
^1  'hn^  Av 

where  Av  is  a  linewidth.  Rado's  gain  is  equivalent  to  Maier 's 
if 

A  =  1. 
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L  =  1. 


and 


e  = 


n 


2 


T  =  l/27rc  Av. 
r 

The  experimentally  observed  gain  for  H2  above  10  atmos- 
pheres  was  stated  as  1.45  x  10  cm/W  for* an  incident  light 
source  of  975  nm.  If  the  gain  coefficient  developed  in  equa¬ 
tion  (AlO)  is  divided  through  by  intensity  then 


2  2 
n 


|3x"l 


Solving  for  1 3x" |  the  result  is 


(A15) 


gcA 

|3x"l  =  - 4  (A16) 

647t-^ 

2 

assuming  n  =1.  The  value  for  3x''  is  the  same  as  that 
used  by  Rado  to  represent  x”* 

The  gain  coefficients  given  in  two  textbooks  more  nearly 
resemble  the  gain  coefficient  in  equation  (A9) .  Pantell  and 
Putoff  [30]  define  g  in  MKS  units  as 


g 


47r  x" 


(A17) 


whereas  equation  (9)  may  be  rearranged  as 


471  |3x"| 


(A18) 


differing  from  Pantell  by  a  factor  of  If  -  1/  they 
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become  the  same.  Yariv  [27]  gives  a  voltage  gain  coefficient 


as 


g' 

=  %  |E 

1^  lx" 

1  9 

(A19) 

The  intensity  gain,  g. 

would  be 

2g'. 

Thus ,  the 

intensity 

gain  coefficient  is 

2g'  =  g  =  -^ 

|E|^  |X”I 

= 

n,c 

|X"I  |E|^. 

(A20) 
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PROGRAM  TO  CALCULATE  THE  LOCAL  PEAK  VALUES  OF  TPE 
RESC^ANT  THIRC-CRDER  NONLINEAR  CIELECTRIC 
SUSCEPTIBILITY, CUt  FOR  A  DIATOMIC  GAS 
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